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Abstract. Spatial and temporal variations in the electron-to-proton mass ratio, /i, and in the 
fine- structure constant, a, are predicted in non-Standard models aimed to explain the nature 
of dark energy. Among them the so-called chameleon-like scalar field models predict strong 
dependence of masses and coupling constants on the local matter density. To explore such 
models we estimated the parameters Ap/p = (ja \, s - H\jb) I H\nb and Aa/a = (a bs - ttiab)/a , iab 
in two essentially different environments, - terrestrial (high density) and interstellar (low 
density), - from radio astronomical observations of cold prestellar molecular cores in the 
disk of the Milky Way. We found that A/i/p = (22±4 stat ±3 sys )xl0" 9 , and |Aa/o-| < l.lxlO" 7 . 
If only a conservative upper limit is considered, then | A/j < 3 X 10" 8 . We also reviewed 
and re-analyzed the available data on the cosmological variation of a obtained from Fe i and 
Fen systems in optical spectra of quasars. We show that statistically significant evidence 
for the changing a at the level of 10~ 6 has not been provided so far. The most stringent 
constraint on lAa/o-l < 2 x 10~ 6 was found from the Fen system at z = 1.15 towards the 
bright quasar HE 0515^-414. The limit of 2 x 10~ 6 corresponds to the utmost accuracy 
which can be reached with available to date optical facilities. 

Key words. Line: profiles - ISM: molecules - Radio lines: ISM - Techniques: radial ve- 
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1. Introduction 

Spatial and temporal variations in the electron- 
to-proton mass ratio, fi = m e /m p , and in 
the fine-structure constant, a = e 2 /(hc), are 
not present in the Standard Model of particle 
physics but they arise quite naturally in grant 
unification theories, multidimensional theories 
and in general when a coupling of light scalar 
fields to baryonic ma t ter is considered dUzan I 
12003b iMartinsI 120081: IChTn et al. I l20o4 . The 



Send offprint requests to: S. A. Levshakov 



light scalar fields are usually attributed to a 
negative pressure substance permeating the en- 
tire v isible Universe and known as dark en- 
ergy (ICaldwellefa/.lll998l) . This substance is 
thought to be responsible for a cosmic acceler- 
ation at low redshifts, z < 1 (Pee bles & Ratal 
120031) . However, scalar fields cause a problem 
since they could violate the equivalence prin- 
ciple which has never been detected in local 
tests dTurvshev et al. 1 120041) . A plausible so- 
lution of this dissension was suggested with 
a so-called 'chameleon' model which assumes 
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that a light scalar field acquires an effective po- 
tential and an effective mass due to its cou- 
pling to matter that stro ngly depends on the 
ambient matter density ( Khourv & Weltman 
2004 iMota & Shaw 112007b lOlive & Pospelov I 
2008). In such a way the scalar field may evade 
local tests of the equivalence principle since 
the range of the scalar-mediated force is too 
short (/leg % 1 mm) to be revealed at the ter- 
restrial matter densities. This is not the case, 
however, for the space based tests where the 
matter density is considerably lower, an effec- 
tive mass of the scalar field is negligible, and 
an effective range for the scalar-mediated force 
is very large (A^s ~ 1 pc). 

Calculations of atomic and molecular spec- 
tra show that different transitions have dif- 
ferent sens itivities to changes in fundamental 
constants (IVarshalovich & Levshakov! 119931 : 
iDzuba et aZ.lll999b iKozlov et alWlQm . Thus, 
measuring relative radial velocities, AV, be- 
tween such transitions one can probe the hy- 
pothetical variation of physical constants. For 
instance, a spatial dependence of fi and a on 
the ambient matter density can be tested lo- 
cally using terrestrial measurements of molec- 
ular transitions and comparing them with radio 
astronomical observations of molecular clouds 
in the disk of the Milky Way (here we have 
a typical difference between the environmen- 
tal gas densities of ~19 orders of magnitude). 
Complementary to these measurements, a tem- 
poral dependence of yu and a on cosmic time 
can be tested through observations of high- 
redshift intervening absorbers seen in quasar 
spectra. 

In this presentation we summarize our re- 
cent results on spatial and temporal variations 
of fi and a obtained at high redshifts z > 
(QSO absorption systems) and at z — (the 
Milky Way disk). 



IChand et all 12004) . This method is based on 
the comparison of wavelengths of different 
transitions in different ions having different 
sensitivity coefficients to the variation of a and 
rising from the same QSO absorption-line sys- 
tem. However, the approach to estimate Aa/a 
on base of different ions seems not to be very 
favorable: it requires many absorption systems 
in order to suppress kinematic effects caused 
by irregular velocity shifts between different 
ions, and its averaging procedure over many 
sight lines and over large range of redshifts 
may smear out a putative weak signal in Aa/a. 
It is clear that in order to obtain Aa/a esti- 
mates at every particular redshift the kinematic 
effects must be suppressed. Such a procedure 
can be realized if lines of only one ion and 
arisin g from the same a tomic levels are uti- 
lized (Levshakov 2004). The corresponding 
method, — the single ion differential a mea- 
surement (SIDAM), — and its applications to 
Fe i and Fe n systems are described below. 

Fei and Fen are represented in optical 
spectra of quasars by several resonance tran- 
sitions having different sensitivity coefficients. 
Under these conditions the accuracy of the 
Aa/a estimates is limited by only two fac- 
tors: (1) uncertainties in the calibration of the 
wavelength scale (see, e.g., Centurion et al., 
these proceedings), and (2) unknown abun- 
dances of iron isotopes in a particular absorp- 
tion system. For quasar spectra taken with 
the VLT, it is possible to reach the accuracy 
in the wavelength ca libration of 30-50 m s _1 
(Molar o et al . 2008a) which translates into the 
error in Aa/a of ~2xl0~ 6 . The boundaries 
on iso topic composition of iron do not excee d 
10" 6 (IKozlov etal.\ 12004 iPorsev et a/.ll2009h . 
Thus, a single absorption system with iron ions 
can provide the accuracy of the Aa/a estimate 
at the level of 10~ 6 . 



2. Cosmological tests (z > 0) 

The cosmological variability of a can be 
probed by different methods. One of them, — 
the many -multiplet method , — was sug- 
gested by Dzuba et al. ( 1999) and extensively 
used in the analysis of quasar metal ab- 
sorbers dWebb et aZ.ll999tlMurphv et a/.l2003l: 



2.1. Fe i afz a b s = 0.45 

Many resonance transitions of neutral iron Fe i 
associated with the absorption-line system at 
Zabs = 0.45 were i dentified in t he sp ectrum of 
HE 0001-2340 bv lD'Odoricol d2007l) . Among 
them there are Fei lines with quite different 
sensitivity coefficients which makes such sys- 
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Fig. 1. Combined absorption-line spectra of Fei as- 
sociated with the Zab S = 0.45 absorption-line system 
towards HE 0001-2340 (normalized intensities are 
shown by dots with lcr error bars). The smooth 
curves show the model profiles. The normalized 
xl = 0.80 (v = 59). The zero radial velocity is fixed 
at z = 0.45207. 

tems suitable for individual Aa/a estimates. 
The sensitivity coeffici ents for Fei transi- 
tions were calculated by Dzuba & Flambaum 
( 2008 ) on our request and revealed to be highly 
sensitive to a change. 

The profiles of the six selected Fei lines 
are shown in Fig.Q] In this case a simple one- 
component model can describe adequately the 
observed intensities. The most probable value 
of Aa/a — (7 ± 7) x 10~ 6 was calculated from 
the Ay 2 -curve presented in Fig. [2] Here we 
combined the Fei AA2967, 3021 A lines hav- 
ing Q * 0.08 and the Fei /U2484, 3441, 3720, 
3861 A lines for which Q 0.03. In Fig. |2] 
the parameter AV is the velocity offset between 
these two groups of Fei lines. In this case the 
estimate of Aa/a is given by 



Aa/a = AV/(2cA<3) , 



(1) 



where AV = Vi — V2 is the difference be- 
tween the radial velocities of these two groups 
of lines, and AQ = Q2 - Qi . 

We expect that with higher spectral resolu- 
tion and higher S/N the accuracy of the Aa/a 
estimate at z a b s = 0.45 can be increased by 
several times since the Fei lines in this sys- 
tem are extremely narrow (the Doppler width 



Av (km . ') 

Fig. 2.x 2 as a function of the velocity difference A V 
between the Fei AA2967, 3021 A lines and ,U2484, 
3441, 3720, 3861 A lines for the one-component 
model shown in Fig. [JJ The lcr confidence level 
is determined by A^ 2 = 1 which gives cta,- = 200 



^ 1 km s 1 ) and have very simple profiles 
(lD'Odoricoll2007l) . 



2.2. Feu afz a b s = 1.84 

The only QSO spectrum which was obtained 
especially with the objective to measure Aa/a 
is the VLT/UV ES spectrum of a bright quasar 
Q 1101-264 dLevshakov et al.\ l2007t) . This 
spectrum was observed with a high resolu- 
tion (FWHM = 3.8 km s" 1 ) and S/N > 100. 
Analyzing the Fen AA26QQ, 2382, and 1608 A 
lines from the absorption system with z a b s = 
1.84, we found AV = -180 + 85 m s~' between 
the Fe n A 1608 A line which has negative sen- 
sitivity coefficient and the Fe n AA26QQ, 2382 
A transitions having almost identical positive 
sensitivity coefficients. In terms of Aa/a this 
shift corresponds to Aa/a = (5.4 + 2.5) x 10~ 6 . 

When new sensitivity coeffic ients for the 
Fe ii lines became available dPorsev et al\ 
2007), we repeated the analysis of this sys- 
tem including the unblended parts of the Fe n 
/12374 A line. The profile of this line is blended 
in the velocity range -30 < V < -10 km s~' 
with telluric absorpti on and because of that i t 
was not considered in Levshakov etal.\ (120071) . 
This line has positive sensitivity coefficient like 
the Fen AA26QQ, 2382 A lines, but its os- 
cillator strength is close to that of the weak 
line Fen A 1608 A. This allows us to con- 
trol possible saturation effects in the strong 
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Fig. 3. Profiles of Fen lines associated with the 
Zabs= 1-84 damped Lyo- system towards Q 1 101— 
264 (normalized intensities are shown by dots with 
lcr error bars). The S/N ratios per pixel (from the 
top panel to the bottom) are 100, 138, 120, 120. 
The smooth curves show the model profiles. Parts 
of the Fen /12374 A line are contaminated in the 
ranges -25 < V - 10 km s~ ] and 4 < V < 13 
km s _I with telluric absorptions and not included in 
the analysis. The positions of the subcomponents of 
the equidispersion 17-component model are marked 
by the dotted vertical lines. The normalized x\ = 
0.87 (v = 300). The zero radial velocity is fixed at 
z= 1.838911. 



lines Fen AA2600, 2382 A caused by unre- 
solved components which may lead to some 
radial velocity shifts between the apparent po- 
sitions of the strong and weak Fen lines. 
The profiles of all four Fen lines (Fig. [3} 
were fitted by a series of mathematical mod- 
els differing both in the number of components 
(from 8 to 17) and in the types of compo- 
nents: either all components have individual b- 
parameters or all components have identical b- 
parameters — the so-called equidis persion de- 
convolution dLevshakov et a/.lll999l) . 

The lcr confidence limit for the veloc- 
ity shift between the Fen A 1608 A and Fen 




Fig. 4.x 2 as a function of the velocity difference A V 
between the Fen A 1608 A and ^2374,2382, 2600 
A lines (shown in Fig. [3} for the 8- and 17- 
component equidispersion models. The lcr confi- 
dence level is determined by A^ 2 = 1 which gives 
<r Av = 90 m s -1 for both represented models. 



X12600, 2383, 2374 A lines was calculated via 
the A^ 2 -curve (Fig.|4]i and by means of Monte- 
Carlo simulations (Fig. 0. In turn, modeling 
of the data points in the Monte-Carlo sim- 
ulations occurred also in two ways: (1) the 
method of bootstrapping residuals, and (2) the 
residuals were modeled by a Markov chain. 
The bootstrap procedure breaks correlations 
between the data points presented in the orig- 
inal QSO spectra (see, e.g., Appendix A in 
Levshakov et al. 2002), whereas the Markov 
chain approximation to the residuals allows us 
to save the original correlation in the simu- 
lated data. In all approaches the negative ve- 
locity shift between the Fe n A 1608 A line and 
the Fen AA26QQ, 2382, 2374 A lines was sta- 
ble reproduced with the most probable value of 
AV = -130 + 90 m s _I (as a lcr limit the most 
conservative value from all trials is taken). In 
terms of Aa/a this corresponds to Aa/a = 
(4.0 ± 2.8) x 10~ 6 . We consider the result ro- 
bust in this respect in disagree ment with what 
claimed by Murphy et al. (2008) who have not 
provided a quantitative analysis of the z = 1.84 
system to support their claim. 



2.3. Feu afz a b s = 1.15 

The most accurate to date estimate of Aa/a 
was obtained for the Fen absorption system 
at z a bs= 115 detected in the spectrum of 
one of the brig htest high-redshi ft quasar HE 
0515-4414 (Reimer sl^ [l9~98). Ouas t et al.\ 
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Fig. 5. Histograms of the simulated samples gen- 
erated by the bootstrapping residuals (panel a) and 
by the Markov chain approximation (panel b) for 
the 17-component equidispersion model. The sam- 
ple size n = 100 in both cases. The smooth curves 
are Gaussians denned by the first two moments of 
the corresponding samples: (4.3 ± 2.0) x 10~ 6 (panel 
a), and (4.5 ± 2.3) x lO" 6 (panel b). 



(2004) used co-added exposures taken with the 
VLT/ UVES and reported Aa/a = (-0.4±1 .9)x 
10~ 6 . IChand et al.\ (12006) re-observed this sys- 
tem with the high resolution and temperature 
stabilized spectrograph HARPS mounted on 
the ESO 3.6 m telescope at the La Silla obser- 
vatory and obtained Aa/a = (0.5 + 2.4) x 10~ 6 . 
We note that in both publications a wavelength 
calibration error was not taken into account, 
and once considered this their errors cr^ a / a 
would be a factor two larger. 

Due to the brightness of HE 0515-4414, 
the S IN ratio in individual exposures of the 
VLT spectra is sufficiently high which allows 
us to pr ocess them separ a tely, i .e. without co- 
adding. iLevshakov et al.\ d2006l) analyzed the 
best quality individual exposures and then av- 
eraged Aa/a values produced from the pairs 
of Fen lines (combinations of the blue line 
Fe n /11608 A and different red Fe n lines). The 
resulting mean value was Aa/a = (-0.07 + 
0.84) x 10~ 6 . However, it was overlooked that 
in this procedure the individual Aa/a val- 
ues from different pairs become correlated. 
Ac counting for t h ese corr elations, as described 
in iMolaro et all (2008b), and using the up- 



dated sensitivity coefficients from lPorsev et all 
(2007) we correct the latter Aa/a value to 
Aa/a = (-0.12 + 1.79) x 10~ 6 . Thus, all Aa/a 
estimates at tubs = 1.15, which are the most ac- 
curate so far, exclude any cosmological varia- 
tions of a at the level of 2 x 10~ 6 . 

3. Local tests (z = 0) 

3.1. Estimate of A/i //j 

Among numerous molecules observed in the 
interstellar medium, ammonia NH3 is of a par- 
ticular interest for the Afi/jj tests due to high 
sensitivity of the inversion frequency (J, K) = 
(1, 1) at 23.7 GHz to a change in fi. Here Afi//i 
= O^obs - MabV/flab- The sensitivity coefficient 
(Q = 4.5 ) of the inversion transition w as calcu- 
lated by dFlambaum & Kozlovll2007l) : 

(Av/v) inv = (y-v)/v = 4.5(A yU /^), (2) 

where v and v are the frequencies correspond- 
ing to the laboratory value of y. and to an al- 
tered ytz in a low-density environment, respec- 
tively. For the rotational frequencies we have 
Q = 1 and 

(Av/v) rot = (v-v)/v = Afi/fx. (3) 

Comparing the apparent radial velocities for 
the NH3 inversion transition, Vi nv , with rota- 
tional transitions, V rot , of other molecules aris- 
ing co-spatially with NH3 one finds 

A/j/h = Q.3(V mt - V im )/c = Q.3AV/C , (4) 

where c is the speed of light. 

The parameter AV in (0 could be consid- 
ered as the sum of two components, AV = 
AV M + AV„, with AVp being the shift due to 
//-variation, and AV n is the so-called Doppler 
noise. The input of the Doppler noise to a pu- 
tative A[ijpL signal can be reduced to some ex- 
tent if the velocity shifts caused by inhomo- 
geneous distribution of molecules and instru- 
mental imperfections are minimized (for de- 
tails, see Levshakov et al. 2009a). 

To probe A///// under different local en- 
vironments, we analyzed at first high resolu- 
tion spectra (FWHM = 25 m s _1 ) of the NH 3 
(J,K) = (1,1) and CCS J N = 2\ - 1 tran- 
sitions obtained with the 100-m Green Bank 
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Tel escope (GBT) by iRosolowskv et al\ d2008l) 
and lRathborne et al.\ (120081) . and moderate res- 
olution spectra (FWHM = 120 - 500 m s 1 ) of 
the NH 3 (J,K) = (1, 1), HC 3 N J = 5 - 4, and 
N2H + 7=1-0 transitions obser ved with the 
45-m N obeyama radio telescope by Sakai et al. 
(2008). These spectra showed a systematic ve- 
locity shift between rotational and inversion 
transitions with the most accurate value of 
A V = 52 + 7 sta t + 13 svs m s" 1 based on th e NH 3 
and CCS lines dLevshakov et al. 1l2008al) . 

Then we carried out our own observa- 
tions of ammonia and other molecules with the 
32-m Medicina, 100-m Effel sberg, and 45-m 
Nobey ama radio telescopes dLevshakov et all 
I2009ah . 

Observations at the Medicina telescope 
were performed with both available digital 
spectrometers ARCOS (ARcetri Correlation 
Spectrometer) and MSpecO (high resolution 
digital spectrometer) which have channel sep- 
arations of 4.88 kHz and 2 kHz, respectively. 
For ARCOS, this resolution corresponds to 62 
m s _1 at the position of the ammonia inversion 
transition (23 GHz) and 80 m s for the rota- 
tional HC 3 N (2-1) line (18 GHz). For MSpecO, 
it is 25 m s _1 and 32 m s _1 at the corresponding 
frequencies. 

The Effelsberg data were obtained with 
the HEMT (High Electron Mobility Transistor) 
dual channel receiver. The measurements were 
obtained in a frequency switching mode, with 
a frequency throw of ~5 MHz. The back- 
end was an FFTS (Fast Fourier Transform 
Spectrometer), operated with its minimum 
bandwidth of 20 MHz and simultaneously pro- 
viding 16384 channels for each polarization. 
The resulting channel widths are 15.4 and 20. 1 
m s _1 for NH 3 and HC 3 N, respectively. 

At Nobeyama we used a low-noise 
HEMT receiver, H22, for the NH 3 observa- 
tions and the two sideband-separating SIS 
(Superconductor- Insulator-Supercondu ctor) 
receiver, T100 dNakaiima et al.\ 120081) . for 
the N2H + observations. Both of them are 
dual polarization receivers. We observed two 
polarization simultaneously. Autocorrelators 
were employed as a backend with bandwidth 
and channel separation of 4 MHz and 4.375 
kHz, respectively. This corresponds to the 



channel widths of 57 m s 1 at 23 GHz, and 14 
ms-'at 93 GHz. 

The most accurate estimate in this case was 
based on the comparison of the radial veloc- 
ities of the NH 3 inversion line with the rota- 
tional transition of HC 3 N 7 = 2-1 observed 
with the 100-m Effelsberg telescope: AV = 
23 ± 4 stat + 3 sys m s -1 . The distribution of ra- 
dial velocity offsets between the rotational and 
inversion transitions is shown in Fig. [6] The 
points represent molecular clouds with sym- 
metric line profiles which can be described by 
a single-component Gaussian model. This re- 
quirement provides a better accuracy of the line 
center measurements. In the selected clouds 
the widths of the emission lines do not exceed 
greatly the Doppler width due to thermal mo- 
tion of material (typical kinetic temperatures 
Tkin ~ 10 K). This ensures that turbulent mo- 
tion does not dominate in the line broadening 
and thus the selected molecular lines sample 
the same kinetic temperature and arise most 
likely co-spatially. 

The difference between the GBT and 
Effelsberg AV values is probably connected 
to uncertainties in the laboratory frequencies. 
The most accurate frequencies for the ammo- 
nia inversion transitions are m easured with the 
uncertainty s v ^ 0.6 m s" 1 dKukolichlll967t 
lHougen|[l972h . The HC 3 N J = 2 - 1 lab- 
oratory frequencies are known with th e er- 
ror s v 3 ms 1 (Mu ller ef a/1l2005h . The 
worst known is the frequency of the CCS 
Jn = 2i — lo transition. The GBT ra- 
dial velocities of CCS were measured with 
the freque ncy v = 22344 033(1) MHz es- 
timated by lYamamoto et al.\ d 19901) from ra- 
dio astronomical observations (1 kHz at 22.3 
GHz corresponds to s v = 13.4 ms _1 ). The 
only laboratory measure ment available give s 
v = 22344.029(4) MHz dLovas et aZ.lll992l) . 
With this frequency we have AV = 6 ± 
7 stat + 54 sys m s _1 . The CD MS catalogue value 
v = 22344.0308(10) MHz dMiiUer et a/.ll2005h 
yields AV = 22 ± 7 stat + 13 sys m s _1 which is 
in line with the Effelsberg measurements based 
on the NH 3 and HC 3 N lines. 

These results show the importance of high 
precision laboratory measurements which are 
expected to provide radio frequencies with un- 
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Fig. 6. Radial velocity offsets, AV, between the 
HC 3 N J = 2 - 1 and NH 3 (J, K) = (1, 1) transitions 
for the data obtained with the 32-m and 100-m tele- 
scopes, and between the N2H + 7=1—0 and NH 3 
(/, K) = (1,1) transitions for the 45-m telescope 
(lo" statistical errors are indicated). /? is the ratio of 
the Doppler 6-parameters: /3 = /?(NH 3 )/fc(HC 3 N), 
or p = £>(NH 3 )/Z?(N2H + ). Given in parentheses are 
ltr errors. The filled circles mark sources with ther- 
mally dominated motions. 



certainties ~ 1 m s 1 to be consistent with con- 
temporary radio astronomical observations. 

If we interpret the most accurate to date 
Effelsberg measurements of AV in terms of the 
electron-to-proton mass ratio variation, then 
A^ = (22±4 stat + 3 sys )xl0- 9 . 

3.2. Constraint on Aa/a 

Complementary to the A^/jU estimate, con- 
straints on a-variation can be obtained from 
the comparison of the atomic fine-structure 
(FS) and molecular rotational transition s 
( Kozlo v et al. 2008; LevshakoveLs3|2008b). 
In this case the velocity difference between the 



rotational and FS lines is sensitive to the com- 
bination of a and /j, F = a 2 /fi: 

AF/F = 2Aa/a - A/j/fi = AV/c , (5) 

where AV = V 10t - Vf s . 

For practical applications we may com- 
pare the FS transitions of carbon [Ci] with 
rotational transitions of 13 CO. The spatial 
distribution of I3 CO is closely traced by 



the 



1997 



si 

Hoi 



C 1] FS lines ( Soaans & van Dishoeck 



0-1 



m s 

cise rest 
J = 1 



llkeda et al\ l2002t iPapadopoulos et al. 
20041) . The rest hyper-fine frequencies of 
low-/ rotational transitions of 13 CO are 
known wi th very high acc uracy, |e v | $ 0.1 
zzoli et al\ 120041) . The most pre- 
frequencies to date for the [Ci] 
transition 49 2160.651(55) MHz 
dYamamoto & Saitolll99lb and J = 2 - I tran- 
sition 809341.97(5) MHz dKlein et al.\ fl998h 
restrict the uncertainties of the line centers to 
e v = 33.5 m s _1 and 18.5 ms _1 , respectively. 
Thus, if we take s v =2 34 m s~' as the most con- 
servative error of the line centering in a pair of 
13 CO-[C 1] lines, then the limiting accuracy for 
AF/F is about 0.1 ppm. 

The analysis of the published results on 
[Ci] and 13 CO low resolution observations 
(FWHM = 0.2 - 1.0 kms" 1 ) of cold molec- 
ular cores gives AV = + 60 sta t ± 34 sys 
m s~' (Levshakov et al. 2009b), which leads to 
a limit |AF/F| < 2.3 X 10" 7 . With the ob- 
tained estimate of A/i/fi, one finds from (0 that 
\Aa/a\ < 1.1 x 10~ 7 . 

This result is already an order of magnitude 
more sensitive than the bound on the cosmo- 
logical a-variation found from absorption sys- 
tems of quasars. The present estimate of Aa/a 
at z = can be considerably improved if higher 
resolution spectra of molecular cores will be 
obtained and the rest frequencies of the [Ci] 
FS transitions will be measured with better ac- 
curacy. 

4. Conclusions 

Our current results of astrophysical tests on 
spatial and temporal variations of fundamental 
constants can be summarized as follows. 

1. The relative radial velocities of the rota- 
tional transitions in HC3N (J — 2-1) 
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and the inversion transition in NH3 (J, K) — 
(1,1) measured with the 100-m Effelsberg 
radio telescope reveal a velocity offset of 
AV = 23 ± 4 stat ± 3 sys m s -1 . This offset 
is regularly reproduced in observations of 
different cold molecular cores with differ- 
ent facilities at the 32-m Medicina and 45- 
m Nobeyama radio telescopes. 

2. Being interpreted in terms of the electron- 
to-proton mass ratio variation, the found 
velocity offset corresponds to Afi/fj. = (22 ± 
4 sta t ± 3 sys ) x 10~ 9 . To cope with negative 
results obtained in laboratory experiments, 
a chameleon-like scalar field is required to 
explain the non-zero Afi value. 

3. A new approach to probe a- variations us- 
ing the fine-structure transitions in atomic 
carbon [Ci] and low-/ rotational transi- 
tions in 13 CO is suggested, and a strong 
constraint on the dependence of a on the 
ambient matter density is deduced at z = 0: 
\Aa/a\ < 1.1 X 1(T 7 . 

4. Another method to measure the cosmolog- 
ical variability of a from resonance tran- 
sitions of neutral iron Fe 1 is discussed and 
applied to the Fe 1 system at z a bs = 0.45. The 
resulting constraint is Aa/a = (7+7)xl0~ 6 . 

5. The Fen system at z a bs = 1-84 is re- 
analyzed and a corrected value of Aa/a is 
obtained: Aa/a = (4.0 ± 2.8) x 10~ 6 . 

6. For the Fen system at z a bs = 115 the 
overlooked correlations between individual 
Aa/a values are now taken into account 
that gives Aa/a = (-0.12 + 1.79) x 10~ 6 . 

7. It is shown that spectral radio observations 
in the Milky Way are at present the most 
sensitive tool to probe the values of funda- 
mental constants at different ambient phys- 
ical conditions. The internal (statistical) er- 
ror of our Effelsberg data is cr AjU ^ = 4 x 
10~ 9 which is about 1000 times lower as 
compared with the error in optical observa- 
tions of quasars. 

8. Since extragalactic molecular clouds have 
gas densities similar to those in the inter- 
stellar medium, the value of A/i//i in high- 
z molecular systems is expected to be at 
the same level as in the Milky Way, i.e. 
Ayu/^-lO -8 — providing no temporal de- 



pendence of the electron-to-proton mass 
ratio is present. 

To conclude, we note that in order to be 
completely confident that the revealed veloc- 
ity shift between rotational and inversion lines 
is not due to kinematic effects in the molec- 
ular clouds but reflects a density-modulated 
variation of A/i//i, new high precision radio- 
astronomical observations are needed for a 
wider range of objects. Besides, such ob- 
servations should include essentially different 
molecules with tunneling transitions sensitive 
to the changes in /j. and molecules with A- 
doublet lines which also exhibit enhanced sen- 
sitivity to variations of yt and a (Kozlov 2009). 
It is also very important to measure the rest 
frequencies of molecular transitions with an 
accuracy of about 1 m s _1 in laboratory ex- 
periments — in some cases the present un- 
certainties in the rest frequencies are larger 
than the errors of radio-astronomical measure- 
ments. In addition, search for variations of the 
fine-structure constant a in the Milky Way disk 
using mid- and far-infrared fine-structure tran- 
sitions in atoms and ions, or search for vari- 
ations of the combination of a 2 //j. using the 
[C n] and/or [C 1] fine-structure transitions and 
low-lying rotational lines of CO would be of 
great importance for cross-checking the re- 
sults. 
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